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Weak Coupling Phase Diagram of the Two Chain Hubbard Model 
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We present a general method for determining the phase diagram of systems of a finite number of 
one dimensional Hubbard-like systems coupled by single-particle hopping with weak interactions. 
The technique is illustrated by detailed calculations for the two-chain Hubbard model, providing the 
first controlled results for arbitrary doping and inter-chain hopping. Of nine possible states which 
could occur in such a spin-1/2 ladder, we find seven at weak coupling. We discuss the conditions 
under which the model can be regarded as a one-dimensional analog of a superconductor. 
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One dimensional (Id) electron systems provide an im- 
portant testing ground for understanding electron corre- 
lation effects. Many methods have been applied to the 
problem of a single Hubbard chain, and there is general 
agreement that the system remains, for repulsive inter- 
actions, in a Luttinger liquid state with gapless spin and 
charge modes [1]. The Id analog of a superconductor, a 
state with one gapless charge mode and dominant pair- 
ing (rather than charge density wave) correlations, does 
not arise in that case. 

Two chain systems are interesting as a first step to- 
wards true 2d materials, and may be relevant for some ex- 
perimental systems [2]. Moreover, on a ladder, statistics 
are more important, since particles can exchange without 
passing through one another. However, the theoretical 
situation in such models is much less clear [3-8]. Re- 
cent simulations suggest that states with dominant pair- 
ing correlations can indeed arise [9]. 

In this letter, we present a systematic weak coupling 
analysis of two Hubbard chains coupled by single par- 
ticle hopping, t±. Unlike previous work, our approach 
is a controlled renormalization group valid for small U 
but for arbitrary inter-chain hopping and filling, n [10]. 
The general methods described here may be applied to 
any system composed of a finite number of Hubbard-like 
chains with weak short range four-fermion interactions. 

The possible phases of such models can be character- 
ized by the number of charge and spin modes which are 
gapless at zero momentum. For an iV-chain system the 
number of gapless charge modes can vary from zero to 
N, and likewise for spin. Remarkably, of the nine possi- 
ble phases for two chains, seven are realized within the 
simple Hubbard model at weak coupling, reflecting the 
proliferation of marginal operators. Denoting a phase 
with x gapless charge modes and y gapless spin modes 
as CxSy, the small U phase diagram for small fixed t± 
is shown in Fig. 1. Particularly noteworthy is the phase 
ClSO, present with purely repulsive interactions (posi- 
tive U). This phase has a spin gap and a single gapless 



charge mode, and is thus the Id analog of either a su- 
perconductor (SC) or charge density wave (CDW). We 
suggest two alternative physical criteria to distinguish 
the two possibilities. Preserving the resistance of the SC 
state to weak impurity scattering requires a slow decay of 
pairing correlations [11]. In particular, if the equal time 
pairing correlation function (A(x)A t (0)) ~ l/|x| K , where 
A = cifC2j., this requires n < k c = 1/3. Alternatively, if 
an array of ladders is weakly coupled via Coulomb inter- 
actions and hopping, the SC state is more resilient with 
k c > 1/2. The traditional requirement [12] of "dominant 
SC correlations" gives k c = 1 (see below). 

Interestingly, the spin-gapped ClSO phase occurs in 
two different regimes (Fig. 1), one for doping, 5 = 1 — n, 
away from half-filling, and the other when the Fermi en- 
ergy coincides with a band edge, fcpi = 0. In the for- 
mer case, pairing correlations develop upon doping the 
spin-gapped Mott insulator at half filling, n = 1, as in 
Anderson's original RVB picture for superconductivity in 
the Cuprates [13]. The critical doping S c at which ClSO 
gives way to a gapless spin state, C2S1 and C2S2, is large 
for weak inter-chain hopping, decreasing from 5 C = 1 for 
small t± to S c = as t± — > 2t. Note that the phase 
C2S2 is the Id analog of a Fermi liquid with all spin and 
charge modes gapless. The presence of the spin-gapped 
state (ClSO) near fcpi = 0, can be attributed to the coin- 
cidence of the Fermi energy with the Van Hove singularity 
at the Id band edge. 
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FIG. 1. Phase diagram at fixed t± < t (not to scale). 



The two chain Hubbard model is described by the 
hamiltonian H = H + H\j , with 
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where c (c^) and <i (e^) are fermion annihilation (cre- 
ation) operators on the first and second chain, respec- 
tively, and a =t, J, is a spin index. The parameters t and 
t± are hopping matrix elements along and between the 
chains, and U is an on-site Hubbard interaction. Eq.l 
has the usual U(l)xSU(2) charge/spin symmetry. 

For weak coupling it is natural to proceed by first diag- 
onalizing the quadratic portion of the hamiltonian. This 
is achieved by canonically transforming to bonding and 
anti-bonding band operators: ipi a = (c a + (—iyd a )/\/2, 
with i = 1,2. In momentum space Hq becomes 



r dp 



(2) 



where ei = t± — 2tcosp and e 2 = —t± — 2tcosp. For 
t± > 2t, the two bands are completely separated. At 
half-filling, the system is then a band insulator, and 
when doped becomes an ordinary spin-1/2 Luttinger liq- 
uid. For t± < 2t, the bands overlap over some range of 
energies. When the Fermi level lies within this region, 
interaction effects must be re-examined in detail. 

It is sufficient to consider the behavior of the sys- 
tem only near the two Fermi momenta kpi, defined by 



= The chemical potential, /i, is fixed by the 
requirement fcpi + kp 2 = nir, where n is the particle 



number per site. The decomposition ip ia « ipRi a e* kFiX + 



(1) fk 



gives, up to a constant, 



H o = ^2 dxvi (ip Ria id x ip Ria - ipl ia id x ip Lia ^ , (3) 



where v t = 2isinfcpi- The allowed 4 fermi interactions 
are highly constrained by symmetry. In addition to U(2) 
invariance, these terms must be preserved by time rever- 
sal, parity, chain interchange, and spatial translation op- 
erations. At generic fillings, the two fermi momenta are 
incommensurate, and the symmetry under translations is 
effectively doubled into independent transformations in 
each band. To delineate the couplings in a physical way, 
we employ the notation of current algebra, 



M iR = -itpRi^tpRn, 



JiR - ^ Ria ^af3fpRif3, 
N Ra (3 = -0Rla"0R2/3, (4) 



where cr denotes Pauli matrices. Left-moving currents 
are defined analogously. There are eight allowed interac- 
tions connecting left and right movers for generic fillings, 
with Hamiltonian densities 



~^Ant — QlpJlR-JlL + 9 2 p J 2 rJ 2 L + 9x P {J\rJ 2 L + ^2R ^1l) + <7lcrJlR - JlL 

t rt- 



2R 
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+<L<t(Jir • J2L + J2R • Jil) + 9t P (L R L L + L R Ll) + (^(Lr • L L + • L[). 



(5) 



Six additional interactions are completely chiral, 



-n 



(2) 



^ip(Jm + Jil) + ^ p (J 2 r + JIl) + ^x P (JirJ 2 r + JilJil) 

+Ai ct (Jir • Jir + Jil - Jil) + A2o-(J2R ■ J2R + J2L ■ J21J + A K(J (Jir • J2R + Jil • J2l)- 



(6) 



The couplings in Eq.6 renormalize "velocities" of various charge and spin modes, and can be neglected to leading 
order in U for what follows. Additional operators are needed to treat umklapp processes at special dopings: 



-?4 3 t = *i(Ml a M 1L + Ml L M 1R ) + u 2 (Ml R M 2h + Ml L M 2R ) + u x (Ml R M 2h + M^M^ + m} r M 1l 
+M 2R Ml L ) + u a (N Ra0 N haP + N Ra0 N[ a(3 ) + u t2 (N Ra(3 N h0a + Nr^N^J. 



(7) 



The single-band umklapp term, Ui, is non-zero only if 
k^i = tt/2. At half-filling the three inter-band umklapp 
terms {u x ,u t i,u t2 ) are non- vanishing. 

The Hubbard model values for the coupling constants, 
obtained from Eq.l, are shown in Table I. To analyze the 
behavior of the weakly interacting system, we employ 
the renormalization group (RG) approach. In the RG, 
short-wavelength modes are progressively eliminated in 
a systematic way, leading to differential equations for 
the renormalized coupling constants which describe the 



physics of the model at longer and longer length scales. 
The flow equations for this system in the absence of u l 
were first obtained in Ref. [14] using conventional dia- 
grammatic methods. The full set of RG equations is more 
directly obtained using current algebra [15]. Away from 
half-filling, they are 

9i P = P(9t P + Y^9ta) ~ au l, 
3 

hp = a(9t P + T^9t*) ~ P u l, 
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9xp — (9t p + jg^toOj 
g 

9la = -ag\ a - -g 2 a + 2(3g tp g ttT , 
g-2a = -PgL ~ 7;9ta + ^g tp g t a, 

9x<T = -gla ~ 2^ta ~ ^9tp9t(T, 

. _ 3 

9tp — 9op9tp + —9o<r9t<T, 

9ta = 9oa9t P + (go P - goa/2 - 2g xa )g ta , 
iii = -2ag lp U!, u 2 = -2/3g 2p u 2 , (8) 

where g { = tt(wi + v 2 )gi, a = (vi + v 2 )/(2v 1 ), (3 = 
(vi + v 2 )/(2v 2 ), g 0p = agip + f3g 2p - 2g xp , and g 0a = 
agi a + /3g 2a — 2g xa . The dots indicate logarithmic deriva- 
tives with respect to the length scale, i.e. g t = dgi/d£, 
where I = In L. 

Eqs.8 are valid until maxjgij} ~ 0(1). To analyze 
them, we employ the following approach. Starting with 
the appropriate initial values (c.f. Table I), we integrate 
the equations numerically. If, as I — > oo, all the couplings 
approach finite values, the procedure is controlled, since 
max{<7j(£ = oo)} becomes arbitrarily small as U — > 0. If 
any coupling diverges, we determine the asymptotic be- 
havior of all the couplings with Eqs.8. As U — > 0, this 
asymptotic behavior is approached arbitrarily closely, as 
the divergence occurs for larger and larger I. While the 
fate of a particular choice of initial {gi} can only be de- 
termined numerically (due to the complexity of the equa- 
tions), the possible asymptotic behaviors when some cou- 
plings diverge can be determined analytically. 

To do so, we make the ansatz <?;(£) = kg l o/(l — kt), 
where 1/k is the scale at which the couplings diverge. 
Eqs.8 then reduce to a set of coupled quadratic equa- 
tions for the {gio}- The search for appropriate solutions 
is considerably aided by the numerical integration of the 
flow equations. After locating a divergence (which fixes 
k), we plot (l — ki)g l versus £, from which g i0 is extracted 
from the intercept with the line I = 1/k. 

Applying this procedure for generic fillings with Hub- 
bard initial values, we found three distinct phases (in 
the regime with both bands partially filled for [7 = 0). 
For a > 4.8, the flows are stable, with fixed point val- 
ues g* a = g* a = g* a = g* p = 0. When 4.3 < a < 4.8, 
the system is singly unstable, with g 2a ,o = — and 
all other g i0 = 0. For more comparable fermi velocities, 
1 < a < 4.3, all the operators except g x(r0 = diverge, 
but in such a way that ag\ a Q = /3g 2a o and gt a o = —^gt P o- 
The behavior for 1/2 < a < 1 is obtained by interchang- 
ing band indices in all quantities. 

The physics of these phases is elucidated through the 
use of abelian bosonization [12,10]. With the convention 
tpR/Lia ex P(«'s/47r0R/Lici)) dual canonical Bose fields 
may be defined as <j> ia = RiQ + (j> Ua and 8 ia = RiQ - 
(f> Lia . They satisfy [<f){x), 6{y)\ = -isgn(x - y)/2. A 



further canonical transformation to (<j),8)i p = [(</>, #);-)- + 
(<Mh]/^2 and {ct>,6) ia = [{</>, 6) n - (<Mh]/^2 yields 
the spin-charge separated Euclidean action 

So = Y,f ^[{dM 2 +{d x 6 iv ) 2 ] + id x e iv d T cl> iv , (9) 

where v = p,a. Using the scheme discussed in the in- 
troduction, the non-interacting system with both bands 
occupied (Eq.9) is classified as C2S2. The large a phase 
found above is also of C2S2 type, though it contains the 
additional (marginal) couplings g ip , g xp , \ ip , \ xp ,\ i(7 , and 
X xa , which makes the behavior highly non-universal [16]. 

In the intermediate state (4.3 < a < 4.8), g 2a becomes 
large and negative. Using bosonization, this interaction 
(neglecting unimportant gradient terms) is 

S 2 «ocg 2lJ [ M 2 cos(V^e 2a ), (10) 

where the coefficient M is cut-off dependent. In the scal- 
ing limit, it is appropriate to expand the cosine and ob- 
tain a true mass M for 8 2a . The resulting phase is there- 
fore C2S1. For a < 4.3, an analogous cosine appears in 
the 1(7 sector, and the asymptotic divergence of gt p and 
g ta is such that the inter-band hopping terms sum to 

Stccgtp / cos(V^</)_ p ) cos(\/27r6li -) cos(v / 27r6'2 t T), (11) 

where <j)± p = (4>\ p ± 4> 2p )/\/2. It is natural to assign 
masses to 6\ a and 8 2a , after which Eq.ll acts as a mass 
for 4>- p . The resulting state has only a single ungapped 
charge mode (with dual fields 8 +p ,(p +p ), and will be la- 
beled C1S0 [10]. It may be characterized by a single 
dimensionless stiffness K +p and a velocity v +p , such that 
the effective action for 8 +p is 

S+p = ^ [ {v +p (d x +p ) 2 +v+l(d T +p ) 2 } , (12) 

Z JX,T 

with a similar dual (K +p — > K+ p ) form for (f> +p . Eq.12 
interpolates smoothly between a CDW (at large K +p ) 
and a superconductor (at small K +p ). The pairing 
exponent k = K +p /2, while the CDW-like correlator 
(n 2 (x)n 2 (0)) c ~ cos[2(A;fi + kF 2 )x]/x 2 / K +" , plus power 
laws (1/x 2 ) at k = 0. We are unable to determine K +p 
and v +p in a controlled way, because they depend on 
the entire crossover from the non-interacting state to the 
C1S0 fixed point. However, heuristic calculations suggest 
that superconducting fluctuations (K+ p ) increase with 
increasing v\ — v 2 and decrease with interactions. In- 
terestingly, the usual power law term at 2k? is missing 
from the density-density correlation function in the CI SO 
phase, as noted by Nagaosa [8] in a similar model, due 
to strong fluctuations of the 9- p and (j> la fields. 
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It remains to discuss the behavior at several "special" 
points in the phase diagram. When kp 2 = tt/2, umk- 
lapp processes imply u 2 7^ 0. For t± > t, this occurs 
with band 1 empty, and one gets the usual C0S1 spin 
density wave. For t± < t, we must consider inter-band 
coupling via Eqs.8. For almost all ratios of the velocities, 
we find that a charge gap develops in band 2, simultane- 
ously suppressing the other potential instabilities, lead- 
ing to a C1S2 phase. Surprisingly, over the narrow range 
0-6 < {3 < 0.85, u 2 renormalizes to zero, yielding instead 
the C1S0 state. The behavior at half-filling is more diffi- 
cult to obtain, because it requires the inclusion of the u x , 
u t i, and u t 2 operators in Eq.7. The RG equations in this 
case are too complicated to reproduce here [15]. Their 
analysis indicates a completely gapped (C0S0) phase, as 
suggested by a large U picture of coupled antiferromag- 
netic Heisenberg chains [15]. 

The final remaining special point occurs when the 
fermi level lies precisely at the bottom of band 1. It 
is outside the scope of conventional RGs, because the 
dispersion in band 1 is quadratic, with the hamiltonian 

H ^ = -I x ^ 4dl+t ^ (13) 

where e = 1 for the quadratic band, but must be taken 
as a small parameter to control the perturbative treat- 
ment. The allowed couplings are g2 P , g2a, Mp, and the 
four interband terms 

= U p (J R2 + Jh2)Jl + W<t(Jr2 + JL2) " Jl 

~\A ~ utiip^ip^ipipipic + h.c), (14) 

where Ji = V'LV'ia and Ji = 4>\ a ^ a p VW 2 - Tne RG 
equations in this case are [15] 

92 P = u 2 /2, g 2rT = -gl - 2u 2 , 

ii p = jm~ € u 2 , u a = -ju 2 a , 

v = ev — 2v 2 — u 2 , 
tit = (e/2 - 2v - 4 7Up - 3 5ct /4 + g p )u t , (15) 

where m = 2m € , 7 = 1/(1 + m~ € ), u± = m e l 2 Ut, 
v = m e v/(4:ir), (g 2 ,u) Pta = (g 2 ,u) Pta /(2irv2) (we have 
taken a momentum cut-off of 1). Because of the relative 
simplicity of Eqs.15, we have been able to analytically 
show an instability for Hubbard initial conditions [15]. 
Analysis of the flow asymptotics indicates that at this 
special point, fluctuation induced attractive interactions 
v < in band 1 populate it with spinless bound pairs. 
Simultaneously, g 2a < creates a spin gap in band 2, and 
u t Josephson couples the two bands, gapping the out of 



phase charge mode to leave a C1S0 phase. This RG anal- 
ysis holds for U < t « 1, but we expect the tendency to 
pairing to increase with e, due to the increased density of 
states for inter-band scattering. Physically, the presence 
of the C1S0 phase at fcpi = can be attributed to the 
Van Hove singularity at the band edge. 

The above results at special fillings are valid only at 
isolated points for infinitesimal U. For finite U, the RG 
suggests that the regions of attraction of these phases 
widen into fans of width St± ~ exp(— ct/U), where c is 
a constant. The non-uniformity of Eqs.15 allows for ad- 
ditional structure within the fan. In particular, because 
the instability is driven by the g 2a term, we expect a 
narrow intermediate wedge of C2S1, as shown in Fig.l. 

It will be interesting to generalize the above calcula- 
tions to three-chain systems (N=3), to help clarify which 
features are particular to even N. 
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TABLE I. Hubbard model coupling constants. The 
g-ology notation is given for comparison with Ref.14. 



Coupling g-ology Hubbard value 

dip Aaaaa/2 - 9aaaa -U/A 

hp 9bbbb/ 2 - 5bbbb -U/A 

9xp <?abab/2 — Sabba —U/A 

dtp Aaabb/2 - flAABB -U/A 

h<* 2p AAAA U 

<?2cr 2g BBBB U 

9x<7 2<? ABAB U 

9t<r 2p AABB U 
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